Breaking arches with vibrations: the role of defects 
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We present experimental results about the stability of arches against external vibrations. Two 
dimensional strings of mutually stabilizing grains are geometrically analyzed and subsequently sub- 
mitted to a periodic forcing at fixed frequency and increasing amplitude. The main factor that de- 
termines the granular arch resistance against vibrations is the maximum angle among those formed 
between any particle of the arch and its two neighbors: the higher the maximum angle is, the 
easier to break the arch. Based in an analysis of the forces, a simple explanation is given for this 
dependence. From this, interesting information can be extracted about the expected magnitudes of 
normal forces and friction coefficients of the particles conforming the arches. 

PACS numbers: 45.70.Mg 



The formation of arches is a common feature when- 
ever large assemblies of solids move collectively. Force 
chains that propagate along a string of particles are able 
to stop the movement, in such a way that the ensemble 
can resist an external pressure as a solid does. Arches, 
defined as arrangements of mutually stabilizing sets of 
particles capable of withstanding external loads [2, 
are the key ingredient for attaining a solid-like structure. 
The likeness of this can happen, for instance, in traffic 
jams [IjHj avalanches of crowds in panic Q and colloidal 
systems {&]. The notion that a common description can 
be given for all those various systems was put forward by 
Cates et al. 0, who called them fragile matter. They 
pointed out that arches can be shattered by the exertion 
of incompatible stresses, meaning forces that act in di- 
rections differing from that of the load supported by the 
arches. 

A large amount of grains is a suitable instance of such 
a situation, amenable to be studied in the laboratory. 
When a dense flow of grains moves through an orifice, 
they are prone to clo g d ue to the formation of an arch 
that blocks the exit [i- 12 1. The properties of these arches 
have been characterized in previous works 13|, Il4j and 
the main findings can be summarized in two points: 1) 
large arches tend to be semicircular in average, i.e. their 
span is twice their height. This implies that there exists 
a direct relationship between all the geometrical features 
(such as the number of beads, the span, and so on). 2) 
There is a considerable number of particles which are sus- 
pended from above the equator, stabilized by frictional 
forces. Concerning this point, in |14[ the angle associated 
to each grain, </>, was defined as the one subtended by the 
two segments connecting the center of the sphere with 
the centers of its two neighbors (see Fig.[T] inset). It was 
reported that 17% of the spheres hang from above the 
equator (</> > 180°) and were called defects. Of course, 
only friction can stabilize a grain in this way. 

The clogging of silos due to arch formation is often 
avoided in practical situations, such as industrial silos, 




FIG. 1: (Color online) Sketch of the experimental set-up. C: 
video camera; A: accelerometer; V: electromagnetic shaker; 
M: motor; S: silo. Inset: A photograph of an arch, showing 
the angle <j> f° r one bead. 



by imparting a vibration to the container [15Hlq . The 
vibration is likely to impose incompatible stresses on the 
arch and to break it. In this paper, we take a closer look 
into this procedure and address the question of whether 
there exists a relationship between the geometric charac- 
teristics of the arches and the external vibration needed 
to break them. This can obviously lead to an improved ef- 
ficiency in the above mentioned industrial processes, and 
also to understand better which are the main ingredients 
affecting the stability of arches. These results could also 
be related with several situations, such as compaction 
dynamics 0,[2^], where the role of arches is crucial. For 
the moment, in this first approach, we have focused on 
the force needed to break an arch by increasing the am- 
plitude of the vibration, at a constant frequency. 

We have set up an experimental device consisting of 
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a two-dimensional symmetrical silo placed on top of a 
magnetic shaker (Fig. [IJ. The silo is made of two sheets 
of transparent polycarbonate (390 mm high and 80 mm 
wide) lined with a conductive coating to avoid electro- 
static charges. The gap between the two sheets is of 1.2 
mm, and we filled it with non-magnetic stainless steel 
beads of 1 mm diameter (in some runs we used brass 
beads of the same size) . At the middle of the silo there 
is a horizontal partition with an orifice of 4.45 mm, di- 
viding the container in two equal compartments. An 
electric motor can rotate the silo half a turn around the 
horizontal axis, through a junction allowing free vertical 
motion of the container. Additionally, a standard video 
camera continuously records the region of the outlet. A 
computer controls all the components with the following 
protocol. Starting from a situation where all the beads 
are in the bottom compartment, the silo is rotated half a 
turn around the horizontal axis, so that the beads start 
to fall through the orifice. The eventual formation of an 
arch that stops the flow is automatically detected by im- 
age analysis, and a photograph of the arch is taken and 
stored. Then, a sinusoidal vibration of 1 kHz frequency 
is switched on, and an amplitude ramp of approximately 
0.09 g/s (g is the acceleration of gravity) is applied to 
the silo. This frequency was chosen because it is an or- 
der of magnitude bigger than the characteristic time that 
it takes for a bead to fall its own diameter from rest un- 
der the action of gravity. Typical amplitudes obtained 
are below the micron range, which is about the size of 
the asperities of the beads. We have checked that the 
residual transversal acceleration is always well below 10 
% of the vertical acceleration. The breaking of the arch 
is detected from the video signal, at which moment the 
maximum acceleration of the sinusoidal forcing T is cal- 
culated from the calibrated input signal. The vibration 
is kept until all the beads are in the bottom half of the 
silo, and the procedure restarts. 

The photographs of the arches (Fig. [TJ inset) can be 
analyzed to obtain </> (the angle associated to each bead). 
The two beads at the end of the arch have not an associ- 
ated angle and indeed, they are not considered to belong 
to the arch. From the particle positions one can also ob- 
tain the geometrical features of the arch (height, span, 
number of beads, and so on). Our first aim, then, was to 
try to establish a relationship between the shape of the 
arches (or any other geometrical characteristic) and the 
force needed to break them, to ascertain whether or not 
some kind of arches withstand the external vibration bet- 
ter than others. Our experiment is tailored to this aim, 
because it provides both the acceleration at the breaking 
point and a photograph of the arch. 

The first noticeable result is that the acceleration 
needed to break an arch, T, decreases with <fi m ax, which 
is the maximum angle found among the beads belong- 
ing to a given arch (Fig. [5]). To further substantiate the 
claim that the grain with <f. 
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FIG. 2: The average maximum acceleration Y imposed by 
the external vibration at the instant of arch breaking, as a 
function of cf> m ax, the maximum angle in the arch. Triangles 
correspond to steel beads, and circles to brass beads. Error 
bars are 95% confidence intervals. 



arch, and that it sets the value of T that the arch can 
withstand, we have inspected two hundred high-speed 
recordings of the breaking process (from the arch forma- 
tion to the arch breaking). We have observed that 64% of 
the arches break at the bead with maximum angle; 12% 
break at a bead touching the one with fimax', another 
12% break at the border (where the value of cj) cannot 
be defined) and the rest break elsewhere. More signifi- 
cantly, we observe that if an arch has a defect, in 95% of 
the cases it breaks just there [2l|]. It is therefore quite 
natural that as the arch breaks at the weakest link, i. e. 
the grain with the largest cf>, then the value of T depends 
on 4>max for each arch. 

Once the dependence of T on the maximum angle found 
in the arch has been revealed, one could examine whether 
r also depends on other properties of the arch (span, 
number of beads, etc.). Recall that all the geometrical 
features of the arches are directly related among them. 
Therefore we can focus on one of them, for instance, the 
number of beads r\ (excluding, as remarked above, the 
two end grains). In Fig. [31(a) we plot T vs. r] for steel 
beads. Clearly, the more beads in the arch, the easier it 
is to break it. But this is just because as rj increases, the 
more likely it is that a larger <f> appears. To illustrate this 
we have calculated the order statistics of <f> m ax , i- e. the 
expected value for the maximum of <fi among a set of size 
r], as calculated from the PDF (probability distribution 
function) of 4> (Fig.[H(c), inset). These values are plotted 
along with the measured ones (see Fig.^b)). As one can 
see, the agreement is good, meaning that the factor by 
which 4> m ax grows as beads are added to the arch can be 
understood just in statistical terms. 

This implies that the angles of the arch are a random 
sample of the PDF of (j>. If the expected value of 4>max 
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FIG. 3: (Color online) (a) Value of F versus the number of 
beads in the arch obtained experimentally (A), and using the 
order statistics as explained in the text (□). (b) Average 
value of <f>max obtained for arches with different number of 
beads r\ (A). Expected value of <f) m ax (□) calculated from 
the PDF of <j>, as a function of the sample size (the number 
of beads rf). (c) V vs. g for small intervals of 4> m ax- A : 
[173°, 175°]; o : [179°, 181°]; o : [187°, 189°]; □ : [191°, 193°]; 
V : [195°, 197°]. Dashed lines correspond to the mean of the 
intervals. Data in this figure correspond to steel beads. 



is used to estimate the average acceleration needed for 
breaking the arch (from Fig. ^ one obtains a prediction 
of r as a function of rj that agrees quite well with the 
experimental results (Fig. [3] (a)). Besides, one can select 
the arches within a small interval of 4> m ax (although this 
reduces drastically the number of samples considered) 
and see whether Y depends on r\ in this subsample. It 
is hard to ascertain any dependence, as can be seen in 
Fig. [3] (c). Hence we can conclude that Y depends on 
the geometrical features of the arch only to the extent 
that these variables enhance or reduce the probability of 
finding a larger <j) max . 

In what follows, we perform an analysis of the forces in 
one bead that can explain the dependence of Y on (j) ma x 
for the case of defects (i.e. beads hanging at <j) > tt). Let 
us consider a bead hanging from above the equator from 
two neighbors in a horizontal, symmetric arrangement 
(see Fig. 0] (a)). Certainly this is not always the case, 
and therefore this explanation can only be deemed as 
approximative. Considering the beads depicted in Fig. 0] 
(a), the normal force N can be shown to be almost equal 
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FIG. 4: (a) Diagram of a defect as considered in the text. 
Forces for the right side of the bead are represented with 
dashed vectors, from left to right: weight, normal force and 
friction, (b) T as a function of 8 (same data that in Fig. [2] 
triangles correspond to steel, and circles to brass) along with 
least-squares fits (dashed lines), (c) Histogram for the angles 
<f> of steel beads (A) and brass (o), in logarithmic scale. The 
dotted vertical line corresponds to (j> c = 190° (steel) and the 
dashed one to <j>c = 192° (brass). In the inset, the whole 
histogram is shown in linear scale. 



(in first approximation) at both sides of the bead, and we 
can work out the force balance for half a bead when the 
arch breaks. At this moment, the friction is mobilized 
and there is an external force due to the vibration F — 
mYg, where m is the mass of one bead. Let us define 
9= ((f)- tt)/2 (see Fig.H(a)). Therefore 



UN cos(0) 



Ar sin (0) + ?M 



mYg 



(1) 



Here [i is the friction coefficient. Assuming that 9 is 
small, which in fact it is, we can write 



r« 9 + — (i-i. (2) 

rag rag 

From this simple model we can infer that there should 
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be a linear relationship between T and 9 (or (j) m ax, be- 
cause it depends linearly on 9). In Fig. [4] (b) we show 
that this linear relationship is fulfilled. From these 
data, corresponding to two different materials, we ob- 
tain N = 6.1 ± 0.5 mg and fx = 0.26 ± 0.08 for steel, 
and N = 4.1 ± 0.5 mg and /z = 0.35 ± 0.08 for brass. 
Unfortunately, it is not easy to provide accurately the 
value of static friction coefficient of the materials used: 
values given usually vary grossly depending on surface 
condition, lubrication, particle shape or sphericity, tem- 
perature, and so on (as we have observed in this experi- 
ment). Widely used values [22| are compatible with the 
obtained ones. In any case, the relevant fact is that the 
proposed argument reproduces the dependence of T on 
and yields reasonable values for [i. 

Even more momentous than the values of the friction 
coefficient is the prediction that the normal force in a 
defect should be about just a few times the weight of one 
bead. One could expect that the average pressure at the 
bottom of a silo, even taking into account the Janssen ef- 
fect, should be much bigger than this. Indeed in a recent 
paper (23j the pressure is proposed to be much bigger 
than the value we give for the normal force in a defect. 
Hence, our result suggests that individual particles form- 
ing a defect are submitted to normal forces smaller than 
the average. It would we worth checking our prediction 
either in experiments (which could be quite involved) or 
in numerical simulations. 

Finally, if we set T = in Eq. ([2]) , we obtain 

<£ c = 7 r + 2/x-^, (3) 

which gives 4> c = 192° for brass and <fi c — 190° for steel. 
These angles correspond to the stability threshold in the 
limit of infinitesimal vibrations. In Fig. |4] (c) the right 
side of the </> PDF is presented. A clear cutoff is obtained 
around the predicted values setting T = in Eq. (|2|). 
This reinforces the validity of the explanation offered. 

We have presented experimental evidence showing that 
arches break at the weakest link, which is the bead in 
the arch which clings to its neighbors with the highest 
angle. Indeed, the value of this maximum angle in the 
arch is the best predictor for the force needed to break it 
when submitted to an external vibration. We provided 
a simple argument to show that a direct relationship is 
expected between the acceleration T and the maximum 
angle (f) m ax, in the case (j) ma ,x > 180°. These results open 
an array of various questions, such as the study of the 
influence of the vibration frequency on the arch stability, 
or the expected smaller time lapse that a big defect would 
endure before the arch breaks. Besides, very hard arches 
should be obtained when defects are absent. 
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